Mode-Locking of mobile discrete breathers 
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We study numerically synchronization phenomena of mobile discrete breathers in dissipative non- 
linear lattices periodically forced. When varying the driving intensity, the breather velocity generi- 
cally locks at rational multiples of the driving frequency. In most cases, the locking plateau coincides 
with the linear stability domain of the resonant mobile breather and the desynchronization occurs 
by regular appearance of type I intermittencies. However, some plateaux also show chaotic mobile 
breathers with locked velocity in the locking region. The addition of a small subharmonic driving 
tames the locked chaotic solution and enhances the stability of resonant mobile breathers. 

PACS numbers: 63.20.Pw,05.45.-a 



I. INTRODUCTION 

Nonlinear lattices provide some of the most interesting 
model systems of macroscopic nonlinear behavior with 
experimental realizations From a theoretical per- 
spective they have been progressively recognized not as 
mere discretizations of nonlinear continuous fields (un- 
avoidable for numerical computations), but as a target 
of interest by themselves, due to the distinctive features 
associated to discreteness. 

Among the variety of behaviors of the lattice nonlin- 
ear dynamics, we focus our attention here on the called 
"intrinsic localized modes" or "discrete breathers" (DB's 
for short). DB's are exact periodic, large amplitude and 
exponentially localized solutions 0. These solutions are 
made possible by the combination of nonlinearity and dis- 
creteness: Nonlinearity allows for solutions out of the lin- 
ear mode bands, due to the frequency dependence of the 
oscillation amplitude. On the other hand, discreteness 
sets an upper cutoff in the band structure and prevents 
the multi-harmonic resonances of DB with extended lin- 
ear modes. These two simple ingredients are enough for 
the existence of DB's, wherefrom the generality and wide 
range of interest of the phenomenon. To visualize an im- 
mobile DB see the upper fig. ^ 

These excitations are not only interesting from a theo- 
retical point of view but with respect to the experimental 
applicability as well, concerning fields as diverse as bio- 
physics (myelinated nerve fibers [lj , biopolymer chains 
3|), nonlinear optics (photonic crystals and waveguides 
4jj), Josephson effect (superconducting devices |(|, 
Bose-Einstein condensates 0) or the physics of glassi- 
ness 0- This makes discrete breathers to be the object 
of a remarkable multidisciplinary interest. 

Unlike localization due to impurities or disorder (An- 
derson), intrinsic localization phenomena support mo- 
bile DB solutions, i.e., exponentially localized oscillations 
where the localization center propagates along the lattice 



as time goes by (see the lower fig. P). Although rigorous 
results apply to the generic existence of non- moving DB, 
much less is known about the conditions for their mobil- 
ity. In this respect, Hamiltonian discrete breathers have 
received (comparatively) more attention than their dis- 
sipative counterparts. However, from the perspective of 
applications to experimental situations, the unavoidable 
coupling (both thermal and non-thermal) of the relevant 
degrees of freedom to a variety of other ones, often de- 
mands to consider open systems where power balance, 
instead of energy conservation, governs the nonlinear dy- 
namics of the lattice. 

In this article we pay attention to the problem of syn- 
chronization and resonant behavior in mobile discrete 
breathers (MB's for short) of the forced and damped 
sine-Gordon lattice (the Id standard Frenkel-Kontorova 
model [UEIE]]), illustrating the effects of time scales 
competition in breather dynamics, in a discrete, dissipa- 
tive and non-integrable context. The two time scales of 
the moving pulse are respectively associated to its fre- 
quency, LOb, and its velocity Vb- Our main results can be 
summarized as follows: 

(i) Locking of the breather mean velocity at some ra- 
tional values of the ratio 2itVb /uib for ranges of pa- 
rameters (coupling, driving strength, ... ) is ob- 
servable. 

(ii) The synchronization of breather velocity is deeply 
rooted in the structural stability of pure resonant 
(to be defined soon) MB's, but it is by no means 
limited to it: The locking island in parameter space 
is generally larger than the linear stability domain 
of the pure resonant breather state. 

(iii) The "extra" locking domain is characterized by 
more complex attracting breathers, sometimes 
chaotic in the breather core, but still keeping a 
locked velocity at large integration times. 
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FIG. 1: Time evolution of two discrete breathers: (a) periodic 
pinned breather, (b) 1/2-resonant (see text) mobile breather. 
The localization center moves a lattice site every two periods 
of the external force 



(iv) When (sub)-harmonic perturbations are added to 
the driving term, the stability of the pure resonant 
MB is enhanced (sometimes substantially), taming 
chaotic dynamics and enlarging the locking step 
size. 

The paper is organized as follows: After this in- 
troductory section we present in section II, in a brief 
but self-contained manner, the relevant and most ba- 
sic aspects of dissipative breather mobility in the forced 
and damped discrete sine-Gordon equation (Frenkel- 
Kontorova model). The definition and characterization 
of (p/g)-resonant mobile DB's along with the extended 
Floquet method for the analysis of their linear stability 
are explained in this section. 

In section III we present our numerical results. We 
show that steps of mode-locking velocity, where {p/q)~ 
resonant solutions exist, are found when varying the driv- 
ing strength. We see that this phenomenon is quite 
general since it is also found for an open set of cou- 
pling parameter values. In section IV we focus atten- 
tion on the unlocking transition, i.e. the transition from 
(p/<7)-resonant locking state to quasiperiodic (irrational 
2nvb/u)b) generic velocity. This transition is character- 
ized as a bifurcation via intermittencies of type I by us- 



ing the Floquet methods reviewed in section II. Section 
V is devoted to the phenomenon of locking enhancement 
by adding small additional subharmonic driving. Finally, 
some concluding remarks are given in section VI. 

II. MOBILE DB IN THE F-K MODEL. 
RESONANT STATES AND THEIR STABILITY. 

Mobile dissipative discrete breathers have been well- 
characterized [T3 in the standard Frenkel-Kontorova 
chain subject to homogeneous periodic driving and vis- 
cous damping (see below). These solutions are attrac- 
tors of the dynamics, and thus they are surrounded in 
phase space by a basin of attraction of initial conditions. 
This fact not only provides fast and accurate numerical 
methods for continuation of generic mobile breathers (in 
contrast with expensive root finding methods for contin- 
uation in the Hamiltonian case), but also guarantees the 
very existence [T^ | of exact moving breathers (in contrast 
also with the Hamiltonian situation, where the stability 
and generality of exact moving solutions is currently an 
issue of debate 0]). 

Generic MB's unavoidably excite extended modes 
(loosely referred to as phonons) which tend to delocal- 
ize energy. However, in dissipative systems, the locally 
excited phonons decay exponentially so that the mobile 
breather keeps on a finite localization length, essentially 
determined by the self-generated phonon dressing. The 
power spectrum (Fourier) analysis of the numerically ex- 
act MB's nicely validates their description as a moving 
source of damped radiation to the extent that the predic- 
tions of the theory exactly match the numerical spectra 
15]. For resonant states (to be defined below) one can 
use Floquet methods in order to perform a more thor- 
ough analysis of these examples of exact non-integrable 
mobility. This will become technically precise along this 
section, after explaining briefly the model and some rel- 
evant issues for the sake of self-containedness. 

The equations of motion of the standard Frenkel- 
Kontorova model subject to a harmonic driving force 
F ac sin(wfct) and a viscous damping a, are, in dimension- 
less form: 

1 • , 

Uj + mij H sin(27rM.,) 

2vr ■>' 

= C (iif+i - 2uj + Uj-i) + F ac sm(ubt) (1) 

where C denotes the coupling (indeed the cou- 
pling/nonlinearity ratio) between neighboring nonlinear 
oscillators Uj(t) of unit mass. 

Two different mechanisms for mobility of DB solutions 
of equation have been observed. 

i) The spontaneous mirror symmetry-breaking of 
pinned discrete breathers (DB), that occurs at 
moderately low couplings, pave the way to mobility 
in a very natural manner, because a moving DB is 
a solution with broken symmetry. Indeed, this sim- 
ple idea is at the origin of a very useful procedure 
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[TM fl6l IT?! to prepare good initial conditions in 
the basin of attraction of exact mobile DB's: As 
described in detail elsewhere 0, adding a small 
perturbation (along the symmetry-breaking eigen- 
vector, often dubbed as depinning mode) to the 
immobile exact DB often evolves asymptotically to 
an attracting moving solution. 

ii) Immobile quasiperiodic DB's have been seen to suf- 
fer from depinning parametric instabilities leading 
to mobile (typically slower) DB's ^| m a range of 
somewhat larger values of the coupling parameter. 
No stable pinned DB coexists in this range with 
moving solutions. 

Once an attracting MB has been precisely determined, 
continuation from it through variations of parameters 
such as coupling or driving intensity, generates a numeri- 
cally continued branch of moving breathers. The contin- 
uation procedure of MB's from an initial one proceeds by 
a very small (adiabatic) change of a parameter, say for 
example A.F ac (or C, or whichever), and numerical inte- 
gration during several periods (Tt) of the external driver. 
The convergence to the breather attractor corresponding 
to F ac + AF ac is guaranteed to be exponentially fast pro- 
vided AF ac is small enough and no bifurcation occurs. 

The range of coupling values where these mobile so- 
lutions are observed is very far from the continuum 
regime; they are roughly in the range from C = 0.5 to 1, 
and its motion has a well defined average velocity. We 
should remark here that these MB's are not related to 
the breathers of the integrable (continuous) sine-Gordon 
model. The continuum limit of Eq. (JIJ is the continu- 
ous sine-Gordon equation under external ac forcing and 
losses, which does not support mobile breather solutions 


A convenient quantitative descriptor of a mobile local- 
ized solution is its translation velocity Vb- In order to 
define precisely this quantity one has to introduce a con- 
tinuous collective variable X(t), naturally interpreted as 
the instantaneous center of localization of energy: 

~ ~rr — s~ (2) 

where e.j = e$ — e^ is the energy density referred to the 
background, i.e.: 

e, = ^ J 2 + ^ I [l-cos(27r Uj )] 

C C 
+ -uj-i) 2 + —(uj+i-Uj) 2 (3) 

and eoo is ej at a site j far away from the exponentially 
localized breather core. The breather velocity is then 
defined as the following long-term average velocity: 

i rto+T 

v b = {X)=\im - Xdt (4) 



Along a single branch of continued MB's, the breather 
velocity defines a continuous but not necessarily smooth 
curve. Let us emphasize that moving solutions of velocity 
Vb posses two characteristic time scales, namely u)b and 
2"KVb- The issue of time scales competition is thus of con- 
cern here, in the sense that one would naturally expect 
the emergence of typically associated phenomena, like 
resonances and synchronization in the behavior of MB's, 
as demonstrated by numerical results to be shown later 
in section III. Now let us introduce some basic notions 
needed for what follows. 



A. Resonant mobile discrete breathers. 

A close view on Eq.Q reveals the basic symmetries of 
the dynamics: 

(a) Lattice translation (homogeneous system), 

C{u 3 {t)} = {u j+1 (t)} (5) 

(b) Discrete time Tb shift (periodic driver), 

T{u J (t)} = {u J (t + Tb)} (6) 

(c) Space-time mirror symmetry (specific to such situ- 
ations) . 

S{ Uj (t)} = {- Uj (t+±T b \ (7) 

This operator combine the mirror symmetry 
lZ(u) = ~u of the local term sin(u) and the center 
symmetry of the sinusoidal driver: S — HT 1 / 2 

Given any solution {u n (t)}, the combined action of 
the symmetry operators C, T and S gives a family of 
solutions generated by the symmetry group. We will sec 
that the existence and structure of this family plays a 
central role in our interpretation of numerical results on 
desynchronization of moving DB's. 

Now we will make precise the notion of resonant DB, 
used in the introductory section. A (p/g)-resonant state 
is a synchronized orbit defined as fixed point of the sym- 
metry element C p T q , i.e. 

u n+p (t + qT b ) = u n (t) (8) 

If p and q have no common divisors, the state is said pure 
resonant. Note that it follows from the definition that 
a (p/g)-resonant MB has a velocity i>& = (uif)/2Tr)(p/q). 
But it is also important to realize that among all conceiv- 
able evolutions {u n (t)} with this velocity, a resonant one 
is certainly very special, for it posses a specific (C p T q ) 
symmetry. 
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B. Linear stability analysis. 



Let us consider a small perturbation {ej(t), £j(t)} of a 



0.04 



given DB solution {uj(t)}. 
tion, we obtain: 



Linearizing around this solu- 



ej + aij + cos{2nu {t))e = C (e j+1 - 2e + e 3 _i) (9) 

An immobile DB of frequency ujf, is a fixed point of the 
operator T. The Floquet matrix maps a basis of the tan- 
gent space ({ej(0), ij(0)}), onto {e,(7b), ej(T&)}, and it is 
given by the Jacobian of T, i.e. D(T). The spectrum of 
eigenvalues of this Floquet matrix gives the linear stabil- 
ity of immobile DB's, allowing the characterization of the 
bifurcations that these solutions experience along contin- 
uation paths, as shown in Refs. 0,^1- In order to use 
these powerful Floquet methods for the analysis of mo- 
bile solutions, they must be periodic and so, one has to 
restrict attention to (p/(7)-resonant MB's. 

Because a (p/g)-resonant MB is a fixed point of the 
operator C p T q , the (extended) Floquet matrix providing 
the linear stability of a (p/<7)-resonant MB is D(C p T q ) = 
VM , where M. is the matrix of the linearized equations 
of motion integrated over q Tf,-periods, and (periodic 
boundary conditions) V is a cyclic permutation matrix 
of p sites. 



Uj(t ) + ejito)} (10) 
{uj(h),ij(to)} + 'PM{e j (t a ),i j (to)} 



The distinctive property of being an attracting solution 
of the nonlinear evolution equation translates into the 
mathematical assertion that an attracting (p/g)-resonant 
state has an associated Jacobian matrix D{C p T q ) with 
(bounded) spectrum inside the complex unit circle: 



sup |/i 



< 1 



(11) 



where /i denotes eigenvalues of the Floquet matrix. 

Eventually, the exit of a Floquet eigenvalue from 
the unit circle signals the destabilization of the [p/q)- 
resonant MB by perturbations along the associated Flo- 
quet eigenvector. In the linear regime these destabilita- 
tions will grow with exponential rate. 



III. THE INTERNAL STRUCTURE OF THE 
LOCKING REGIME 

First of all, we briefly review the method used to gen- 
erate mobile discrete breathers (MB). We begin by gen- 
crating an immobile discrete breather, starting from the 
anti-continuum limit (C = 0). Initially, we have used the 
same parameters as in Il5| : F ac = 0.02, uii, = 0.1 x 2tt 
and a = 0.02. As we increase C adiabatically, the dis- 
crete breather solution remains as an attractor of the 
dynamics |15|. enabling us to find breathers at different 
values of C, by continuation. 
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FIG. 2: Velocity of the MB vs. F ac for C = 0.55. A very 
narrow plateau 3 appears. The rest of fixed parameters are 
u b = 0.1 x 2tt and a = 0.02 



We continue the breathers until the first pitchfork bi- 
furcation [Tol Hil rTi| , which connects one-site breathers 
with two-site breathers via asymmetric ones. The lo- 
calized eigenmode, responsible for this instability is 
asymmetric, and can be used to "depin" the static 
breather. Then, we generate a MB by perturbing the 
static breathers along this mode with an amplitude /x. 
As in jl5l Il6| , we found MB if the perturbation is larger 
than some critical value fi c . Unlike Hamiltonian systems, 
the velocity reached by the MB is independent of [i for 
[i > (jl c . This method allows us to produce MB in a 
wide range of the coupling parameter C. Two kinds of 
MB's have been observed depending on C: for C in the 
interval [0.5,0.89] only induced MB's exist in coexistence 
with static breathers. However, in the C range [0.89,0.97] 
spontaneous MB's appear as attractor solutions, coexist- 
ing whith the induced MB's. Static breathers are not 
found in this range. Hence, we have decided to focus our 
research on some points in those regions. The selected 
values are C — 0.55 and C = 0.75, in the first region and 
C = 0.94 in the second. 

In order to check the dependence of these MB's on 
the driving force, we fix C and then we vary F ac . 
Some general features emerge in all cases. The sim- 
ulations show the appearance of steps with velocities 
t>5 = (uj\ 3 /2tt)(p/ q). Recall that these velocities are the 
velocities of mode- locked MB. But this does not guaran- 
tee that the MB inside the locking step is (p/q)-resonant, 
so we must check the periodicity of these MB's inside 
these steps. Another empirical observation is that the 
limit value of F ac before the destruction of the MB, in- 
creases with C. 

The results with C — 0.55 are sketched in fig. [21 This 
figure shows the full F ac range in which MB's with a 
definite velocity exist. One can see an extremely nar- 
row step at velocity Vb = (wb/27r)(l/3). The Poincare 
section of these solutions reveals that they are pure reso- 
nant MB's. Moreover, the linear stability analysis reveals 
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FIG. 3: Velocity of the MB vs. F ac showing two mode-locking 
plateaux (g and for C = 0.75, the values of 7 and cut are 
the same that in fig. [5] 



that they are linearly stable inside the whole step inter- 
val, i.e. [0.01568,0.01581]^. This type of solution, i.e. 
(l/3)-resonant, has been found in the whole range of C 
values that we investigated. 

For C = 0.75, the results of the velocity-force curve are 
represented in the Fig. |3J In this case we find a bigger 
step at velocity v b = (w 6 /2tt)(1/2) ([0.03177, 0.042] Fac ), 
and, again, a very little one at Vf, = (a>b/27r)(l/3) 
([0.01534, 0.01538] Fa J. In the whole interval of this last 
step, the MB is (1/3) -resonant and linearly stable. On 
the other hand, the bigger step has a more complex 
structure: In the range [0.03177, 0.03831] Fac , the MB 
is periodic with period 2Tb and linearly stable. These 
(l/2)-resonant solutions can be continued in coupling 
parameter from C ~ 0.71 up to C ~ 0.84. However 
at F ac ~ 0.03831 the MB suffers a transition of period 
tripling (a Floquet eigenvalue and its complex conju- 
gate cross the unit circle at angles 27r/3 and — 2tt/3), 
and rapidly goes to a chaotic state via a period doubling 
cascade, for F ac ~ 0.03833. The MB remains in this 
chaotic state with commensurate velocity, and around 
F ac ~ 0.042, the solution leaves the step. 

Finally, at C = 0.94 we obtain the curves represented 
in the Fig. There we can observe a very rich behavior 
of the breather velocity: different steps at velocity values 
of (ct> 6 /27r)(l/3),(ct> 6 /27r)(4/9) and (w b /2vr)(2/3), as well 
as an evident hysteresis. The latter is somehow typical of 
underdamped systems; it implies the coexistence of (at 
least) two different MB attractors, with the same model 
parameters which has been reported previously |l5l Il9| . 
The dynamics in these steps is simple: periodic and lin- 
early stable, with no bifurcations inside the plateau. Also 
a smaller plateau with velocity (ct>&/27r)(l/2) is found, al- 
though in this case only quasiperiodic MB's exist. For the 
sake of completeness we mention that the upper branch 
can be continued to lower values of the coupling and 
connects with MB's generated by the depinning mode, 
whereas the lower branch belongs to the spontaneous 
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FIG. 4: Same as figures |5] and [J] now for C = 0.94. Several 
mode-locking steps are visible. Note, also the hysteresis in 
the curve. 



MB. 

Summarizing, the step with velocity (w£,/27r)(l/3) is 
found in all the range of coupling where MB's seem 
to exist. The range for the step solutions with vi, = 
(u>b/2n)(l/2) starts at C = 0.71 since below this coupling 
a MB with regular motion does not exist, in the parame- 
ter range studied. At higher C, we observe other locking 
steps like Vb — (w{,/27r)(4/9) and Vb = (uJb/2n)(2/3) for 
C = 0.94. Although we have not been exhaustive vary- 
ing in the parameter C we can certainly find these steps 
(and others) in the whole range of coupling parameters, 
where DB are found. 

Some of the observed locking plateaus coincide with 
the stability interval of the corresponding (p/g)-resonant 
MB, but in other cases (notably the 1/2 locking step 
shown in Fig. the resonant state destabilizes inside the 
plateau and the new attracting solution (with locking ve- 
locity Vb) is more complex: either periodic with a larger 
period, or even chaotic, as revealed by the computed 
(largest) Lyapunov exponent. These types of complex 
behaviors inside a locking plateau are known to happen 
for a single driven-damped anharmonic oscillator [2(| , as 
well as for moving discommensurations |2l|. In this re- 
gard, the result for discommensurations can be reducible 
to the single-particle case using collective variable ap- 
proaches. In our case, the breather internal degrees of 
freedom, together with the X (breather center) variable, 
do not allow for a straightforward reduction to single- 
particle behavior. 



IV. UNLOCKING TRANSITION 

A typical route for the transition from a periodic state 
to a quasi-periodic or chaotic one in dissipative systems 
is that mediated by intermittencies 22]. Intermitten- 
cies occur whenever the behavior of a system seems to 
switch between qualitatively different unstable periodic 
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orbits or behaviors (periodic, quasi-periodic or chaotic) 
even though all the control parameters remain constant 
and no external noise is present j^. Depending on the 
type of Floquet instability of the periodic orbit respon- 
sible for the bifurcation (crossing the unit circle at +1, 
at two complex conjugate eigenvalues or at —1) intermit- 
tencies are classified into intermittencies type I, II, or III, 
respectively |22| . 

In this section we study the unlocking transition of 
the (p/q)-resonant MB and characterize it as a transi- 
tion from periodicity to quasi-periodicity driven by the 
regular (periodic) appearance of type I intermittencies. 
This mechanism for unlocking transitions was already ob- 
served in the purely dissipative dynamics of ac driven 
modulated structures of the FK model [24j ■ Intermitten- 
cies of type I are also known to be responsible for the 
depinning transition of discrete solitons (discommensu- 
rations) of the underdamped FK model |21j . 

From now on in this section we will focus our attention 
on the unlocking transition at the left edge of the 1/2 
locking step of Fig. 01 (C = 0.75) which occurs at F ac = 
0.03177. The extended Floquet analysis of the (1/2)- 
resonant MB close to the edge reveals that an eigenvalue 
of the Floquet matrix D(£T 2 ) approaches the value +1 
from the interior of the complex unit circle along the real 
axis, and leaves the unit circle at the transition point. 
The eigenvector associated to this Floquet eigenvalue is 
exponentially localized at the breather center. Out the 
step, the (l/2)-resonant MB is thus linearly unstable. 

In order to visualize the effect of this instability we plot 
in Fig. [3 the breather center (out of the step but very 
close to its edge) X and its velocity X in a reference frame 
moving with the locking velocity (a>(,/27r)(l/2). One can 
clearly see in figureElthat the breather center remains for 
some time intervals in laminar regimes (of locked veloc- 
ity), interrupted by sudden jumps of very short duration. 

We have numerically checked that these laminar 
regimes correspond to the discrete family of equivalent 
unstable continuations of the (l/2)-resonant MB, related 
one to each other by symmetry operations C, S and T. 
Therefore, the destabilizing Floquet eigenvector pushes 
the weakly unstable resonant MB u towards its (equiva- 
lent) closest member of the family (that turns out to be 
5~ 1- u), which is also unstable, and so on. The duration 
of each laminar regime, which diverges at the bifurca- 
tion point, is proportional to (/i — 1) , where /i is the 
unstable Floquet eigenvalue. 

The computation of the Power Spectrum of X, i.e, 



S(u,) = 



X(t)e iuJt dt 



(12) 



for the attracting MB out of the locking plateau (see fig- 
ure E}, reveals the new frequency uji n t = fi — 1 associated 
to the intermittencies, and further confirms that they ap- 
pear at regular time intervals, so that the attracting MB 
out of the step is quasiperiodic. 

This scenario of unlocking mechanism is confirmed to 
happen for other numerically obtained plateaus of mode- 
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FIG. 5: Upper panel: X — Vi oc kt, i.e. X respect to a moving 
frame that moves with Vb = (ub/2n)l/2 just below the left 
edge of the i step of the Fig[3] Jumps correspond to the 
intermittencies described in the text. Lower panel: X in the 
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FIG. 6: Power spectrum S(u>) of X showing clearly a 
quasiperiodic behavior. 



locked velocity. It appears that the desynchronization by 
regular intermittent phase shifts pulses is quite a generic 
phenomenon. 



V. THE EFFECTS OF SUBHARMONIC 
PERTURBATIONS OF THE DRIVING FORCE 

In the previous sections, we have shown how com- 
plex the dynamical response of the system can be. 
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This includes, among others, resonant MB's which are 
destabilized via intermittencies when the parameters are 
changed. Inside the resonant step, the particular behav- 
ior of the oscillators forming the breather can be complex 
(highorder periodic or chaotic) . However it does not pre- 
vent the breather to have a definite mean velocity, com- 
mensurate with the external frequency. 

Our goal in this section is to enlarge the regions of 
the parameter space for which those resonant steps ex- 
ist. The structure of the quasiperiodic attractor in the 
vicinity of the resonant steps gave us some indications 
about the procedure to follow. We must reinforce the 
laminar phase (correspondingly inhibit the appearance 
of intermittencies) applying a perturbation subharmoni- 
cally related with the original ac force: 

F(t) = F ac sm(oJbt) + Asin (— t + A) (13) 

n 

where A is small compared to F ac , n is a positive integer 
and A represents a phase shift between both terms. 

Such kind of perturbation has been proved to be effi- 
cient to stabilize linearly unstable periodic orbits. It has 
been used in systems with few degrees of freedom l25ll2fl 
as well as in solitons j^l and even in experiments |28j|. 
In our case we achieve to lower the onset of the resonant 
steps significantly. In particular, we focus our attention 
on the step 1/2 in C = 0.75 which presents the richest 
phenomenology, and choose n — 2 in i|13|) and A = 
(other A values have been used and the results obtained 
are essentially similar). The main results are summarized 
in the phase diagram of Fig. For example the onset 
of the resonant step is reduced from F ac — 0.031777 to 
F ac ~ 0.013 when A increases from to A ~ 0.001. 

This effect adds to another that takes place inside 
the step and that is related with the control of chaos. 
A chaotic attractor like the one developed by a period 
doubling cascade inside the 1/2 resonant step for values 
F ac k, 0.03833 (with A = 0) is formed by a dense set of 
unstable periodic orbits of different periodicities but all 
of the same velocity. The addition of a suitable pertur- 
bation is able to stabilize one of these unstable periodic 
orbits. The taming of chaotic states in nonantonomous 
dynamical systems by the addition of harmonic pertur- 
bations is a quite general phenomenon |25|, 123, |23 > an d 
we observe it in the case of these chaotic breathers with 
locked velocity. 

To quantify this behavior we compute the largest Lya- 
punov exponent [3(j of the MB solution at a fixed F ac = 
0.04 as a function as the perturbation strength A. We 
start at a chaotic mode- locking MB. As soon as we ap- 
ply the perturbation, a significant decrease of the Largest 
Lyapunov exponent is observed until a narrow region of 
quasiperiodicity is reached (Largest Lyapunov exponent 
= 0). Then, a sequence of periodic and quasiperiodic 
solutions follows and finally, a broad region of periodic 
solution, with the perturbation period (see Fig. [BJ. Note 
that all this is attained with a A two orders of magnitude 
smaller than F ac . 




0.013 0.031777 0.04 F 



FIG. 7: Schematic phase diagram of the behavior with a sub- 
harmonic perturbation of strength A vs. unperturbed driving 
force amplitude F ac . Shaded region refers to more complex 
behavior (chaotic, quasiperiodic or higher resonances than 
1/2) but also mode- locking. The arrow stands for the path 
followed in Fig. |H] 



0,008 




-0,1?'' 



4e-05 8e-05 0,00012 

X 

FIG. 8: (upper) Largest Lyapunov exponent (£) for F ac = 
0.04. (lower) Poincare section of X at intervals 2Tb, so a 
single point is indicative of a periodic solution in the mode- 
locking step at its corresponding A value. 



At fixed C, one can tune, by varying A, the desired 
periodic state. This method is extremely robust against 
changes in C. The pure resonant MB solution can be ex- 
tended for a wide range of the coupling C. Using an ap- 
propriated A we have extended the 1/2-resonant solution 
from [0.71,0.84] to [0.5,0.84] in the coupling parameter. 
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VI. CONCLUDING REMARKS AND 
SUMMARY 

We have studied mobile discrete breathers in the un- 
derdamped Frenkel-Kontorova model. Periodic (or p/q- 
resonant) solutions have been found for a large range 
of parameter values (as C and F ac ). These states are 
structurally stable as they are a consequence of the syn- 
chronization between the two time scales of the mobile 
breather. This synchronization (or its absence) results in 
a very rich dynamics of the mobile breather solutions, in- 
cluding quasiperiodic and chaotic ones. Localized chaotic 
behavior has been previously observed in static breathers 
|3l| but, here we show that it is compatible with mode- 
locking motion of the breather center (core). 

One important issue in the study of localized discrete 
excitations is the possibility of their dynamical descrip- 
tion by reduction to a system with few degrees of freedom 
(or collective coordinates). Up to our knowledge, two 
collective coordinate scheme have been developed with 
mobile breathers. The first approach described in [3^ 
has got very fruitful results for discrete [2l| and con- 
tinuum |33| solitons, as well as for continuum breathers 
[34l l35| . This approach uses as starting point for the 
calculations the continuum model solution and fails far 
from this limit. 

The second approach [36( is the geometric counterpart 
of the first. It works with mobile breathers in a hamil- 
tonian framework, which prevents its direct use in the 
dissipative case considered here. However, the possibil- 
ity of its use in non-hamiltonian contexts, under certain 
technical conditions, is still open as was briefly pointed 
out in ref. .37] 

Although analytical methods for a collective coordi- 
nate approach to mobile dissipative discrete breathers 



have not been yet developed, our numerical results 
strongly suggest that such an approach could give an ac- 
curate account of most of the observed phenomenology. 
We hope that the results presented here could encourage 
for the search of those methods. 

Finally, we remark the important role that a second 
harmonic in the forcing plays in the dynamics of the 
MB. We observe that the presence of a small subhar- 
monic driving term enhances (p/g)-resonant solutions. 
Increasing the subharmonic amplitude i|13|) from zero, an 
initially chaotic MB can be tamed, above some thresh- 
old, obtaining a non-chaotic one (quasiperiodic or (p/q)- 
resonant). The feasibility of the experimental implemen- 
tation of different wave forms for driving forces could 
falicitate the observation of these solutions since they are 
more robust and stable. 

In summary, we have shown that mode-locking motion 
of breathers are stable solutions for an important exam- 
ple of non-linear lattice such as the Frenkel-Kontorova 
model. We have characterized this mode-locking solu- 
tion as well as the mechanisms for the unlocking. Finally 
we have applied a simple method to extend this kind of 
solutions to a large range of values in parameter space. 
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